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ABSTRACT 


A  new  set  of  region  corrections  has  been  generated 
for  LASA  which  more  adequately  covers  the  seismically 
active  areas  of  the  earth  than  did  previous  sets.  The 
new  set  contains  183  calibration  nodes  versus  105  nodes 
on  the  set  used  in  the  SQAC/LASA  system  throughout  1972 
and  the  first  half  of  1973.  Each  node  contains  the 
average  relative  travel- time  anomalies  for  each  of  Ihe 
21  subarrays  at  LASA  which  are  valid  over  an  area  sur¬ 
rounding  the  location  of  the  node.  The  corrections 
were  generated  from  more  than  1800  events  using  as  a 
primary  source  of  data  the  time  shifts  resulting  from 
the  crosscorrelation  process  in  the  SDAC/LASA  Event 
Processor. 


TABLE  OF  CONTENTS 


Page  No. 

ABSTRACT 

INTRODUCTION  x 

The  Current  System  2 

The  New  Node  and  Anomaly  Set  5 

CONCLUSIONS  9 

ACKNOWLEDGEMENT  10 

REFERENCES 


LIST  OF  FIGURES 

Figure  No.  Figure  Title  Page  No 

1  Inverse  Velocity  Space  Deployment  4 

of  Region  Correction  Nodes  Sup¬ 
plied  by  IBM 

2  New  Inverse  Velocity  Space  Deploy-  7 

ment  of  Region  Correction  Nodes 


LIST  OF  TABLES 


Table  No. 

Table  Title 

Page  No. 

IA 

New  Region  Correction  Nodes 

11 

IB 

New  Travel-Time  Anomalies  in 
Hundredths  of  Seconds 

19 

-  vii- 


INTRODUCTION 


In  array  processing,  beamforming  is  the  principal 
method  of  detecting  teleseismic  signals.  Basically, 
the  method  requires  that  suitable  relative  time  delays, 
appropriate  to  the  array  elements,  be  determined  for 
the  location  at  which  the  beam  is  aimed.  To  calculate 
the  delays  some  earth  model,  or  travel  time,  and  distance 
relationship,  must  be  assumed.  However,  since  the 
earth  is  not  known  well  enough  to  predict  relative 
travel  times  to  within  0.1-0. 2  seconds,  an  additional 
procedure  must  be  followed.  When  a  teleseismic  signal 
is  recorded  by  the  various  array  elements,  the  relative 
arrival  times  can  be  observed.  By  comparing  the  ob¬ 
served  relative  times  to  the  theoretical  relative  times 
computed  from  a  known  epicentral  location  and  from  the 
adopted  model,  travel- time  anomalies  can  be  calculated 
in  the  manner  suggested  by  Chiburis  and  Dean  (1965, 
p.  2).  The  beam  delays  can  then  be  corrected  to  com¬ 
pensate  for  the  anomalies  such  that  acceptable  signal 
alignment  of  future  events  is  achieved. 

Our  experience  has  shown  that,  although  the  array 
anomalies  are  reproducible  (that  is,  constant)  for  events 
from  any  particular  geographic  region,  they  vary  signif¬ 
icantly  from  one  region  to  another  (Chiburis  and  Dean, 
1955;  Chiburis,  1966;  and  Chiburis,  1968).  This  fact 
requires  that  a  set  of  anomalies  be  calculated  for 
-very  region  where  an  array  beam  is  aimed. 
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The  data  processing  system  for  the  Large  Aperture 
Seismic  Array  (LASA)  at  the  Seismic  Data  Analysis  Center 
(SDAC),  as  programmed  by  IBM  under  Contracts  F-19628-C- 
0198  and  F-19628-68-C-0400 ,  employs  105  regions  (or  nodes, 
to  emphasize  their  point  nature)  for  which  corrections 
are  assumed  known.  The  LASA  Detection  Processor  (DP) 
has  300  beams  spatially  distributed  so  that  the  signal 
loss  between  beam  centers  is  generally  less  than  3db. 

Each  set  of  beam  delays  has  been  modified  by  corrections 
which  were  determined  by  a  tripartite  interpolation 
scheme  using  the  105  nodes.  Clearly,  a  good  correction 
set  which  offers  adequate  coverage  is  critical  to  all 
of  the  array  processes  of  beamforming  and  location. 

This  report  discusses  the  current  set  of  nodes, 
the  interpolation  scheme,  and  the  establishment  of  a 
new  set  containing  nearly  200  nodes,  each  consisting  of 
21  subarray  corrections  as  determined  from  more  than 
1800  well-recorded  events. 
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The  Current  System 

The  Event  Processor  (EP)  of  the  SDAC/ LAS A  system 
has  as  one  of  its  processing  options  a  correlation 
scheme  whereby  each  sub  array  time  series  is  crosscorre- 
lated,  with  lags,  with  the  array  beam  pointed  at  the 
DP  location.  When  the  signals  are  large  enough  for 
this  process  to  work,  a  set  of  time  shifts  (lags)  is 
determined  from  which  a  final  bean  is  formed.  The 
correlated  time  shifts  are  modified  by  the  assumed- 
known  corrections,  after  which  a  velocity  and  an  azimuth 
are  determined  and  a  corresponding  location  is  computed. 
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Should  the  crosscorrelation  process  fail  because 
of  poor  signal  power,  of  signal  decorrelation,  or  of 
insufficient  lags  for  a  badly  misaligned  detection  beam, 
the  EP  invokes  a  beampacking  process.  In  this  process, 
a  set  of  18  beams  is  iteratively  formed  around  the 
detection  beam  and  the  beam  having  the  largest  power  is 
selected.  Each  of  the  beams  in  the  pack  has  predicted 
delays  modified  by  the  interpolated  corrections. 

Between  January  1971  and  May  1972,  the  SDAC/LASA 
system  operated  with  the  crosscorrelation  method.  As 
shown  by  Ahner  (1973),  beampacking  is  far  faster  and 
yields  results  as  good  as,  and  often  better,  than 
crosscorrelation.  Therefore,  since  May  1972,  beam¬ 
packing  has  been  the  principal  method  of  the  EP. 

The  current  node  set  and  the  tripartite  tie  points 
generates  by  IBM  are  shown  plotted  in  U-space  in  Figure 
1.  For  events  occurring  between  nodes  and  within  a 
triangle,  a  proportional  3-way  linear  interpolation  is 
made  to  obtain  the  subarray  corrections.  It  can  be 
seen  that  about  1  dozen  of  the  105  nodes  cover  the 
region  where  U  <  0.04  sec/km  (beyond  the  P-range  dis¬ 
tance)  .  Interpolating  between  c^ "rections  for  P  and 
(for  example)  PKP  is  invalid;  in  this  case,  it  would 
probably  be  better  to  apply  no  time  corrections  at  all 
to  events  beyond  the  P-range.  The  remaining  nodes  in 
Figure  1  provide  dense  enough  coverage  only  in  the  upper 
left  quadrant  (to  the  northwest  of  LASA,  from  southern 
Alaska  through  the  Aleutians,  Kamchatka,  The  Kuriles,  and 
Japan,  to  the  Marianas).  About  one  fourth  of  the  105 
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nodes  inadequately  cover  the  remaining  three  quadrants. 
As  an  example  of  the  poor  coverage,  consider  the  Azores 
Islands  Region  and  the  N^-th  Atlantic  Ridge  in  the 
northeast  quadrant  (azimuth  65-85°)  at  U-radius  of 
about  0.06  sec/km.  As  seen  in  Figure  1,  all  corrections 
for  this  area  would  be  interpolated  between  a  node  loca¬ 
ted  in  Mcrroco  and  two  nodes  in  Eastern  Montana.  It  is 
not  surprising  that  events  occurring  in  this  region 
reported  by  SDAC/LASA  often  display  misalignments 
greater  than  0.5  seconds  with  the  result  that  the  array 
beam  signal  is  3-6db  smaller  than  it  should  be.  More 
importantly }  events  from  this  area  which  are  near  but 
above  an  achievable  detection  threshold  of  LASA  are 
probably  being  missed  entirely  in  DP  because  of  improp¬ 
erly  corrected  delays.  Referring  again  to  Figure  1, 
there  are  only  four  nodes  covering  the  Asian  continent, 
and  no  node  at  all  for  -he  Hawaiian  region.  Beams 
directed  in  these  areas  are  almost  certainly  being 
poorly  steered. 

This  lack  of  coverage  in  certain  regions  and  improp¬ 
er  coverage  in  others  prompted  us  to  undertake  a  re¬ 
design  of  nodal  placement  and  a  redetermination  of  a 
correction  set  as  completely  as  possible. 

The  New  Node  and  Correction  Set 

In  an  earlier  report,  Chiburis  and  Dean  (1965) 
showed  that  a  comparison  of  array  time  shifts  deter¬ 
mined  by  visual  observation  and  by  cross  correlation  of 
very  good  signals  could  yield  results  within  0.1-0. 2 


seconds  of  each  other.  Therefore,  the  time  shifts 
resulting  from  the  crosscorrelation  process  in  EP 
between  February  1971  and  May  1972  were  used  as  a  pri¬ 
mary  source  of  data  from  which  to  calculate  region 
corrections . 

The  procedure  used  was  to  punch  cards  from  the 
Event  Tape  containing  the  epicentral  parameters  and 
the  21  relative  time  shifts.  The  events  were  associated 
with  those  reported  by  the  Earthquake  Research  Lab  (ERL) 
of  the  National  Oceanic  and  Atmospheric  Administration 
(NOAA) .  From  the  ERL  hypocenters  and  the  crosscorrela- 
ted  time  shifts,  travel- time  anomalies  were  computed 
using  the  Herrin  (1968)  travel-time  tables,  which  are 
the  tables  used  in  the  SDAC/LASA  system.  More  than 
40,000  anomalies  were  computed  in  this  way  and  regional¬ 
ized.  Because  the  crosscorrelation  process  occa¬ 
sionally  misaligns  a  subarray  signal  by  a  full  cycle, 
and  because  it  often  misaligns  a  signal  by  0.2-0. 3 
seconds,  a  considerable  number  of  the  time  shifts 
must  be  visually  read  from  the  plots  of  the  event  in 
order  to  establish  good  corrections. 

The  regions  into  which  the  events  were  placed 
define  the  183  nodes  shown  in  Figure  2.  The  tripartite 
tie  points  were  determined  with  a  knowledge  of  the 
seismicity  between  nodes  and  with  a  view  toward  weighting 
seismically  active  regions  more  than  relatively  aseismic 
regions  for  the  three-point  interpolation  scheme.  As 
• een  in  Figure  2,  the  quadrants  are  about  equally 
covered  with  nodes.  It  is  estimated  that  nearly  80% 
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Figure  2.  New  Inverse  Velocity  Space  Deployment  of 
Region  Correction  Nodes. 


of  the  highl/  seismic  regions  between  30°  and  100°  fm' 
LASa  are  adequately  covered  by  the  183  nodes. 

The  implementation  of  the  new  region  correction 
set  is  expected  to  increase  the  performance  and  effi¬ 
ciency  of  the  SDAC/LASA  system  as  regards  detection 
threshold,  location  accuracy,  EOC  analysis,  and  the 
reprocessing  of  misaligned  signa  >. 

Table  IA  lists  the  Node  Number,  Region  Name, 
Geographic  and  U-Space  location,  the  Number  of  Events 
which  were  used  to  determine  the  corrections  for  that 
node,  and  the  LASA  velocity  and  azimuth;  Table  IB  lists 
the  set  of  anomalies  (in  hundredths  of  seconds)  for 
the  nodes.  For  implementation  as  corrections,  the 
values  in  Table  IB  must  be  reversed  in  sign. 

Additional  events ,  occurring  in  new  regions  or  in 
regions  where  only  a  few  events  were  available,  are 
being  analyzed  on  a  time  available  basis  and  will  pro¬ 
vide  for  periodic  updates  and  improvements  to  the  new 
region  correction  set. 

A  report  on  the  performance  of  the  new  set  relative 
to  the  original  and  current  set  will  be  issued  as  soon 
as  enough  data  have  been  accumulated. 
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CONCLUSIONS 


Our  experience  with  operating  the  SDAC/LASA  system 
led  us  to  believe  that  a  large  percentage  of  signals 
were  not  being  as  well  aligned  as  possible  in  the  beam¬ 
forming  process.  This  misalignment  would  cause  the 
system  to  miss  signals  from  events  above  the  achievable 
detection  threshold  as  well  as  seriously  mislocate  the 
detected  events.  A  preliminary  analysis  of  the  current 
placement  of  105  nodes  and  of  the  region  correction  set 
m  the  system  confirmed  that  many  of  the  seismic  regions 
between  30°  and  100°  were  not  properly  calibrated  or 
were  not  calibrated  at  all. 

A  new  set  consisting  of  183  nodes  was  developed. 

A  region  correction  set  for  the  21  subarrays  at  LASA 
was  generated  from  more  than  1800  events,  using  as  the 
source  of  data  the  time  shifts  resulting  from  the 
crosscorrelation  process  in  EP. 

The  implementation  of  the  new  region  correction 
set  is  expected  to  improve  the  detection  threshold,  to 
yield  more  accurate  locations,  and  increase  the  effi¬ 
ciency  of  the  SDAC/LASA  system. 
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